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Abstract. An elegant formula for coordinates of carbon atoms in a unit cell of a 
single-walled nanotube (SWNT) is presented and the potential of neon (Ne) inside an 
infinitely long SWNT is analytically derived out under the condition of the Lennard- 
Jones potential between Ne and carbon atoms. Specific heats of dilute Ne inside 
long (20, 20) SWNT are calculated at different temperatures. It is found that Ne 
exhibits 3-dimensional (3D) gas behavior at high temperature but behaves as 2D gas 
at low temperature. Especially, at ultra low temperature, Ne inside (20, 20) nanotubes 
behaves as lattice gas. A coarse method to determine the characteristic temperature 
T c for low density gas in a potential is put forward. If T ^> T Cl we just need to use the 
classical statistical mechanics without solving the Shrodinger equation to consider the 
thermal behavior of gas in the potential. But if T ~ T c , we must solve the Shrodinger 
equation. For Ne in (20,20) nanotube, we obtain T c » 60 K. 
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1. Introduction 

Since the discovery of carbon nanotubes [I,, the peculiar electronic and mechanical 
properties of these structures have attracted much attention [3 El El- Experiments 
have also revealed that they can also be used to store hydrogen jH] and other gases [Hj. 
Many physicists expected that gases in nanotubes or nanotube bundles may display 
novel 1-dimensional (ID) behavior as a consequence of the remarkable aspect ratio of 
the length of tubes to their radius. The group led by Cole [Zj and other researchers jH] 
have theoretically studied properties of gases in nanotubes or nanotube bundles. One 
of the most fantastic properties they found is specific heat of dilute gas inside single- 
walled carbon nanotubes (SWNT's) as a functions of temperature: With temperature 
increasing it shows the thermal behavior changing from ID, 2D, and 3D. However, there 
are still two question arisen: Do ID and 2D behaviors always exist for dilute gas inside 
SWNT only if the temperature is low enough? Is classical statistical mechanics (CSM) 
sufficient to deal with this problem? 

In Ref. jH] , we has given a brief answer to the first question. Here we will exhibit the 
full calculations and continue to discuss the second question. In recent papers, Siber et 
al treat the specific heat of dilute He atoms adsorbed in interstitial channels and grooves 
of carbon nanotube bundles ^H] ■ Here we partially follow their aprroach. 

This paper is organized as follows: In SecEJ we give the expressions of the 
coordinates of carbon atoms in a unit cell of a SWNT. In SecQ we analytically calculate 
the potential of Ne in a carbon nanotube. For example, we calculate the potential of Ne 
in (20,20) tube. In SecEl we list the related formula of thermodynamics in statistical 
mechanics to calculate the specific heat. In Sec|Hl we calculate the specific heats of 
Ne inside (20,20) at different temperatures without considering 9, z effects (9, z are 
defined in SecHJ). In Sec® we calculate the specific heats of Ne inside (20,20) at low 
temperature. In Sec0 we discuss why we do not use CSM to calculate the specific heat. 
In this section, we present a coarse method to estimate the characteristic temperature, 
below which the CSM can not be used. In SecJHl we discuss the reliability of our results 
and give brief conclusions. 

2. the coordinates of carbon atoms in a unit cell of a SWNT 

A SWNT without two caps can be constructed by wrapping up a single sheet of graphite 
such that two equivalent sites of hexagonal lattice coincide [TT]. To describe the SWNT, 
some characteristic vectors require introducing. As shown in Fig^ the chiral vector C^, 
which defines the relative location of two sites, is specified by a pair of integers (ni,^) 
which is called the index of the SWNT and relates to two unit vectors ai and a2 
of graphite (C^ = nia.\ + ri 2 a 2 ). The chiral angle #0 defines the angle between a x 
and Cfc. For (nun?) nanotube, 6q = arccos[ — 2ni+n2 =]. The translational vector 

T corresponds to the first lattice point of 2D graphitic sheet through which the line 
normal to the chiral vector passes. The unit cell of the SWNT is the rectangle 
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denned by vectors and T, while vectors &i and a 2 define the area of the unit cell 
of 2D graphite. The number N of hexagons per unit cell of SWNT is obtained as a 
function of ri\ and n 2 as N = 2(nf + n\ + n\n-i) j 'du, where dn is the greatest common 
divisor of (2n 2 + ni) and {2n\ + n 2 ). There are 2N carbon atoms in each unit cell of 
SWNT because every hexagon contains two atoms. To denote the 2N atoms, we use 
a symmetry vector R to generate coordinates of carbon atoms in the nanotube and is 
defined as the site vector having the smallest component in the direction of CV From 
a geometric standpoint, vector R consists of a rotation around the nanotube axis by 
an angle \1/ = 2n/N combined with a translation r in the direction of T; therefore, R 
can be denoted by R = 1 t) . Using the symmetry vector R, we can divide the 2N 
carbon atoms in the unit cell of SWNT into two classes: one includes N atoms whose 
site vectors satisfy 

A, = /R- [/R- T/T 2 ]T (Z = 0,l,2,---,iV-l), (1) 
another includes the remainder N atoms whose site vectors satisfy 

B, = IR + B - [(/R + B ) • T/T 2 ]T 

-[(/R + B )-C,/C 2 ]C, (Z = 0,l,---,iV-l), (2) 

where B = (^ol r o) = ( y^~p — ~ a cc cos(#o — f ) j represents one of the nearest 
neighbor atoms to A . 

If we introduce a cylindrical coordinate system (r, 9, z) whose z-axis is the tube 
axis. Its r#-plane is perpendicular to z-axis and contains atom A in the nanotube. r 
the distance from some point to z-axis, and 9 the angle rotating around z-axis from an 
axis which is vertical to z-axis and passes through atom A in the tube to the point. In 
this coordinate system, we can express Eqs.((TJ) and (0) as 

A l = {p,W 1 lT-[lr/T]T) (/ = 0,l,2,---,iV-l), (3) 



and 



B, = (p, M + * - 27r[**±*°], It + r - 

Z7T 1 



(/ = 0,1,2,...,JV-1), (4) 

where p = In Eqs.(^)-(0J), the symbol [• • •] denotes the largest integer smaller than 
e.g., [5.3] = 5. 



3. the potential of ne inside carbon nanotubes 

To obtain the potential of Ne inside the nanotube, we firstly consider another simple 
system shown in Fig. |2J Many atoms distribute regularly in a line form an infinite atom 
chain and an atom Q is out of the chain. The interval between neighbor atoms in the 
chain is T, and the site of atom Q relative to atom can be represented by numbers c\ 
and c 2 . We take the Lennard- Jones potential U(Rj) = 4e[(a/Rj) 12 — (a/Rj) 6 ] between 
atom Q and atom j in the chain, where Rj is the distance between Q and atom j, and 
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e = ^e c e ne , o = (cr c + cr ne )/2 with e ne = 35.6 K, a ne = 2.75 A, e c = 28 K and a c = 3.4 
A[7J^|. We calculate the potential between atom Q and the chain as 

U QC = 4e[ ( x 12 [/ 6 (c 1 , c 2 ) - a G U 3 ( Cl , c 2 )], (5) 

oo 

where £4(ci,c 2 ) = J2 u , i-, 2 , ztts = 1, 2, • • •) which can be calculated through 

[yc\+nl ) +c 2 J 

n=— oo 

the following recursion [T3] : 



7rsinh(27TC2/T) 



^l(Cl, C 2 ) C2 T[cosh(27rc 2 /T)-cos(27rci/T)] : 

C/fe+i(ci,c 2 ) = -l/(2kc 2 )dU k /dc 2 . 



(6) 



The (20, 20) tube, for example, with infinite length can be regarded as 2N = 80 
chains. Thus the potential of any point Q inside the tube can be calculated as 

2N 

U(r, 0,z) = J2 Uqc, (7) 

i=l 

where (r, 9, z) is coordinates of Q in the cylindrical coordinate system. As an 
approximation, we neglect the potential varying with z and 9 because we find that 
it is much smaller than the potential varying with r through calculations, and fit the 
potential with U(r) = 4e[(^) 10 — (-^) 4 ], where p = 13.56 A is the radius of the tube, 
£ = 390 K, and a = 2.63 A (see also FigOJ). Moreover, we simplify it as 

U(r ) = i -(£)*], r<p-a, 

oo, r > p — a. 

4. The basic formulae of thermodynamics in statistical mechanics 

When we deal with the thermal behavior of dilute Ne in nanotube, we will begin with the 
free energy T = —Tin Z, where T is temperature and Z = e~ En / T = tr exp(— Tt/T) 
is the partition function. We have let the Boltzmann factor be 1 [Tl]. E n is the 
eigenvalues of the Schrodinger equation: Hip n = E n ip n . From the free energy, we 
can derive out the specific heat c v = —T(d 2 J r /dT 2 ) v = ( En '~( En > ; where {E 2 ) = 
J2 n E 2 n e- E ^/Z a nd(E n ) = j: n E n e- E ^/Z. 

We call above statistical mechanics ordinary statistics (OS) and list the 
corresponding formula in CSM: The free energy JF C/ = —Tin Z c i, where Z c i = 
J e~ E( - p ' q ^ r dT is the partition function and here the prime means we integrate only 
over the regions of phase space which correspond to physically different states of 
particles. The specific heat c vc i = —T(d 2 J r d/dT 2 ) v = ^ E , where (E 2 ) = 
j> E 2 e- E ^/ r dT/Z cl and (E) = J' Ee~ E( - p ' q ^ r dT / Z d . 

Under the sufficient high temperature, T converge to the T c \- In other words, the 
applicable domain of T is larger than that of T c \- Therefore, the conclusion derived from 
T is much more reliable than that derived from T c \. We will use OS in the following 
discussions if we do not make special statement. 
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5. Specific Heats of Ne gas in an approximate potential 

Because we consider the dilute Ne, we can neglect the interaction between Ne atoms 
and write the single particle Schrodinger equation J5] as Hip = Eip, where H = 
-^V 2 + U(r) and ip = (Pe^ m9+KZ \ It follows that 

E = ^ + E m (/ceR,m = 0,±l,±2,...), 

H{r)cj> = E m <f>, (9) 

where m is angular quantum number and E m is the corresponding energy. Setting 
r = (p — cr)£, e = an d r] = a/p, Eqs.(jHJ) and Q are transformed into 




!0 _ ( u ) 4 1 t < 1 



(10) 



and 



ii) 



£Q ( d 2 

(l-r,)^d£2 ^ fdf g 

If we let \<f) = ^2 n a n \Xn), we will obtain the secular equation 

det(H jn - E m S jn ) = 0, (12) 

where H jn = J XjT^iOXn^, and S jn = J* XjXn£d£- H we let x« = J\ m \(v n £), where 
Ji m | (£) is the m-th order Bessel function of the first class and v n is the n-th zero point 
of Bessel function [IS], we can calculated E mn , (m = 0, ±1, ±2, • • • ; n — 1, 2, 3, . . .) from 
Eq.(H2]). 

If there are TV Ne atoms inside the tube, we have the free energy T = —NTlnZi, 

ft 2 ft 2 

where Z\ = J2 mn e~ Em "^ T e~ 2 ' tT dn. We can easily obtain the specific heat per 
atom is 

_ T3 2 T _ 1 ~ (-Enm) 2 n on 

" ~ MdT 2 ~ 2 + T 2 ' 1 J 

where (£ nm ) = £ ^^wr and (^ m ) = S ^Vwr • 

In FigJ3J the symbols "V" reflect c„ varying with the temperature T, which implies 
Ne atoms inside (20,20) tube behave as 3D gas at high temperature (specific heat 
approaches to 3/2) and 2D gas at low temperature (specific heat is 1). Therefore, we 
can naturally assume that all atoms in the valley of potential U(r) at low temperature, 
i.e. lie on the shell S* with radius g = p[l — (5/2) ly/6 77] . 

At high temperature, the Ne atoms maybe evaporate from the nanotube, but our 
model cannot recover this effect because our potential is infinite near the surface of 
carbon nanotute. 
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6. Specific Heats at low temperature 

Now we consider the thermal property of Ne inside the nanotube at low temperature in 
detail. We assume that all atoms lie on the shell S* with radius g — p[l — (5/2) 1//6 ?7] 
at temperature lower than some critical temperature T/ c which depends on the energy 
difference between the tow lowest energy of radial excitations. From Eqs.(JSJ)-((Z|) we can 
easily calculate the potential U s (v,z) on S*, where v = g6. The Hamiltonian of single 
particle can be expressed as H' = — f^(J^ + -§zj) + U s (v, z). 

In fact, U s has periodic structure. If we denote Qi = (27ig/N,r), et 2 = (0, T) and 
ji = hoci + /2CK2, we have U s (r + 7;) = U s (r), where h, I2 G Z and r = (v, z). On the 
one hand, we have the Bloch's theorem [T7] : 

which suggests that $(/«, r) = V. a(K + Gj)e^ K+G ^' r , where Gj = jxfli +32^2 with 
ji,j 2 G Z, /3i = (N/g,0) and f3 2 = (-tN/(Tq), 2tt/T). From Eq. pty . we obtain the 
secular equation 

det(Hij - EJ l:i ) = 0, (15) 

where Hij = g(/e + Gj) 2 5ij + U lh U tj = ± j no e i( - G '- G ^ r U s {r)dr and tt = |a x x a 2 |. 

On the other hand, Periodic boundary condition along the circumference of the 
shell S* suggests that we just need to consider the first Brillouin zone which consists of 
k = (m v / g, k z ) where m v G Z, < m v < N and k z G K, < k z < 2n/T. 

From Eq. (jl5|) . we can calculate the energy E mvtKz , and then the free energy 
T = —AfTlnZi [T%] . where Z\ = J2m v e~ Emv ' K ^ T dn z . Moreover, the specific 
heat per atom is 



MOT 2 T 2 ' 1 ' 

N—l N—l 

where (E>) = E !? ,T E mv ^ e - E ^/ r dn z /Z and (£' 2 > = E E^ e^^dK 

m v =0 m v =0 

In FigJHl the symbols "V" reflect c v varying with the temperature T, which implies 
Ne atoms inside (20,20) tube behave as the lattice gas [T^j at ultra low temperature 
(specific heat is 0) and 2D gas at low temperature (specific heat approaches 1). There 
is no ID gas inside (20,20) tube, which is quite different from our usual notion. 



7. why not use CSM? 

At the sufficient high temperature, T converge to the T c \- If a carbon nanotube, for 
example (20, 20) tube, just a geometric tube with diameter d ~ 2 nm, the characteristic 
temperature T c = h 2 /(2irmd 2 ) = 0.04 K. We just need to use the CSM if T >> %. 
However, here the carbon nanotube provides potential to Ne. In the potential, it is 
difficult matter to calculate the characteristic temperature is not a simple matter. For 
example, consider some atoms in a harmonic potential U{r) = ^muj 2 r 2 . We can easily 
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calculate the specific heat per atom c v = (3u; 2 /T 2 ) exp(— u/T)/[l — exp(— uj/T)] 2 [14j. 
Here we have set H — 1. We know c v = 3 if T >> uj and this is the classical case. Thus 
we can set T c = u (In the whole form T c = hLo/ks)- Above discussion suggests that T c 
depends on the potential. We propose a coarse method to obtain T c which is the least 
root of the following equation: 

U m + T C = U(R 1 ) = U (R 2 ) 

R 2 -R 1 = h/V^m% ' [ ' 

The physical meanings of parameters in Eq. (jl7|) are shown in the FigEl 

For the potential U(r) = \muj 1 r 1 in above example, U m = 0, Ri = —R2, R2 — R\ = 
2i? 2 , using Eq. ffTTj) we arrive at T c = ^-uj ~ uj (the whole form is T c = ~ hw/ks)- 
For (20,20) nanotube, the approximate potential is U{r) = 4e[(-^) w — (^:) 4 ], with 
p = 13.56 A, e — 390 K, and a = 2.63 A. We can numerical solve Eq.((T2j) and give 
T c ?z 60 K. If T >> T c (e.g. T > 3T C = 180 K), the classical case is expected and 
CSM is valid. But CSM can not be used if T ~ T c (e.g. T < 2T C = 120 K), especially 
T < T c = 60 K. 

In FigHJwe give results obtained from OS and CSM. The information shown in FigQ] 
agrees with our coarse estimate to T c very well. Therefore, using the CMS is insufficient. 
It is necessary to consider the quantum mechanics. In other words, quantum mechanics 
must be used when studying the behaviors of gases in the nanometer scale. 

8. discussion and conclusion 

We notice that the potential U is a function of (r, 9, z) and U(r, 9 + 2nji/N, z + j\T + 
j 2 T) = U(r,9,z), where N,t,T are parameters of a carbon nanotube, and ji, j'2 are 
two integers. Under this symmetry, we must use a generalized Bloch's theorem (see 
appendix) to find the eigenvalues of Schrodinger equation. We can also obtain the 
secular equation of energies. Although we have 16 Pentium IV computers, but the data 
and computing time exceed our computers' ability if we calculate the specific heat with 
temperature varying from to 400 K. But, fortunately, we find that the variation of 
potential with 9, z is less than 30 K, much smaller than that with r. In SecEJ we neglect 
the effect of the variation of potential with 9, z. Indeed, this effect can be neglected under 
high temperature in terms of intuition. At high temperature, the values of c v approach 
to that calculated from CSM. We show them in Fig|7[ in which the triangles are values 
neglected the effect of 9, z and the line is the result calculated from CSM considering 
the effect of 9, z. It suggests our intuition is reliable. 

We know the effect of 9, z is important under the low temperature. In SecEl we 
consider this effect in detail. Because the shape of U is very sharp and deep with r, the 
particles are confined in the bottom of U(r) at low temperature. Thus we can consider 
the effect of 9, z but neglect the effect of r at low temperature. 

Therefore, we believe that we have grasped the main factors under at high and low 
temperature and our qualitative conclusions are reliable. 
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The interactions (include the contribution from the polarization of the pi-electrons) 
between Ne and C atoms are intermolecular forces, van der Waals type interactions. As 
the theoretical study, the simplest form of van der Waals force is the Lenard- Jones 
potential. We believe that the eventual results based on the Lenard- Jones potential 
are qualitatively reasonable. Maybe, many researchers Ej use it because of the same 
reasons. Indeed, the more detailed study needs to consider the Gay-Berne potential [2U] 
or other complex potentials. We are considering to overcome this difficult problem by 
using numerical simulations. 

It is necessary to notice that we just consider the dilute Ne and the interactions 
between Ne atoms are naturally neglected. Therefore we do not describe the phase 
transition of Ne in SWNT's. If we intend to deal with the phase transition, we must 
include the interactions between Ne atoms and use the methods of quantum field theory, 
such as the renormalization group approach [21], etc. 

In conclusion, we calculate the specific heat of Ne in (20,20) SWNT at low and 
high temperature and find the dimensional crossover of thermal behavior in this system. 
Especially, the dilute Ne gas exhibits as 2D behavior at low temperature or 2D lattice 
behavior at ultra low temperature. The simple physical picture is that all atoms are 
confined on a shell blow a critical temperature. 
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Appendix: Generalized Bloch's theorem 

In SecEl we calculate the specific heat with neglecting the effect of 9, z\ and in Sec® in 
fact, we neglect the effect of r. In other words, we do not obtain the specific heats in 
a consistent way from low to high temperature. Indeed, the potential U is a function 
of (r, 9, z) and satisfies U(r,9 + 2irji/N, z + j\T + j 2 T) = U (r, 9, z), where ji, j 2 are two 
integers. If a potential satisfies this condition, it is called the spiral symmetric. Under 
this spiral symmetry, we must use a generalized Bloch's theorem to find the eigenvalues 
of Schrodinger equation and then calculate the specific heat in a consistent way [2*2"] . 

Above all, we go over traditional translational symmetry and Bloch's theorem |23j . 

The hamiltonian of a system with the translational symmetry is expressed as 






(19) 
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J(Rj)J(Ri) = J(R t + Rj-) = J{Rj + R,) = ^(R^R,), (20) 

and 

JtfllWW = -^J(Rj)V 2 f(r) + J(Rj)V(r)f(r) 
= W(Ri)/(r) + V(r + R,.)/(r + R,) 

= h^V 2 + 7(r)]/(r + Rj) = HJ(K 3 )f(r). (21) 

Therefore {j7"(Rj), if} is the set of conserved quantities. In this case, an eigenfunction 
of the Hamiltonian must be an eigenfunction of the translational operators, i.e., 
i7"(Rj)^(r) = ^(r + Rj) = \(Rj)ip(r) if Hip(r) = Eip(r). Furthermore, the electron 
density must be periodic, i.e., \ip{v + Rj)| 2 = |VK r )| 2 - It follows that 

|A(R,)| 2 = 1. (22) 

But from Eq. (|20|) . we know 

A(Rj)A(Rj) = A(Rj + R). (23) 

The solution of Eq. (J2H|) under the constraint Eq. (|^2*j) is A(Rj) = e lK ' Rj . Thus we have 
Bloch's theorem: 

V>(r + Rj) = e iK ' R ^(r). (24) 

Next, we will turn to the spiral symmetry and generalize the Bloch's theorem. The 
hamiltonian with the spiral symmetry in cylindrical coordinates is expressed: 

h 2 d 2 d d 2 d 2 

2fi^dr 2 rdr r 2 d0 2 dz 2 ^ ' ( ) 

V{r, 6, z) = V{r, 6 + jrf, z + j x r + j 2 T), (26) 

where = 2ir/N, Nt = MT, N,M G N, ji,j 2 e Z. Define vectors r = (9,z), 
Rj — (ji^,ji T + 32T) in the space [0, 2ti) x K, and operators J~(Rj) which act on 
a function /(r;r) as J r (R i )/(r;r) = /(r;r + Rj). It follows that J(Rj)J(Ri) = 
J(Ri + Rj) = JiRj+Ri) = J(Ri)J(Rj) and J(Rj)V(r; r) = l/(r; r) J(Rj). 
Otherwise, 

^ R iX^ + r) = (^ + ^)^)/(r; r). (27) 

Thus 

J(Ri)H = HJ(R j ). (28) 

We can obtain generalized Bloch's theorem analogizing the method to obtain the 
traditional Bloch's theorem: 

ip(r; r + R,-) = e^ R ^(r; r), (29) 
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where k = (I, k), I = 0, 1, • • • , N - 1, < k < 2n/T. 

If we set Ql = (At), ck 2 = (0,T), ft = (AT 0), /3 2 = (—tN/T,2-k/T), 
Gj = ji/3i +32^2 = (jiN — j2Nr/T, 2irj2/T), then we can expand the wave function as 
the superposition of planar waves: 

^(r; r) = ^ C>(r)e i(K+G J> r . (30) 
j 

In cylindrical coordinates, Eqs.(j2HJ) and (j3T?j) are expressed 

1>(r, e + ^,z + j x r + j 2 T) = e^^nr+nT)]^ ^ ^ (31) 

Vfc(r, M) = £ C7 i ^^(r)e«^-*^^* + **^*I. (32) 



JU2 



If Expanding (fi K (r) with normalized orthogonal basis {XnZre^)) ^ £ N} we obtain 
the secular equation from Eq.(j22J) and Shrodinger equation 

det[(l~C mim 2j 1 j2nn' -^nlK^jimi^jama^nn'}] ^ 5 (33) 
where 7~~Cj 1 j 2 rnirn2nn' 1jij2nn'$jimi$j2m2 ^jij2Tnim2nn' > Ijijznn' 2u ./o [,^2 "I - 



2/i JO 

f 2 *df) C T 'dz f p rr/rv-,JrV7rr.fl. 7 rW-„..rrV«[0i-^i)-^ 2 ^]^+ <2 " (J V m2)z 



^ - P**"-^ 3 - (« + ^) 2 ] Xn , K (r)rdr, and 



Ujij 2 m im2 nn' = £f f£* d6 J Q dz f£ rdrXni K (r)U (r, 6 , z)xn'i K (r)e l 
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Figure 1. The unrolled honeycomb lattice of a SWNT. By rolling up the sheet 
along the chiral vector C^, that is, such that the point A coincides with the point 
corresponding to vector Ch, a nanotube is formed. The vectors ai and &2 are the real 
space unit vectors of the hexagonal lattice. The translational vector T is perpendicular 
to C/j and runs in the direction of the tube axis. The vector R is the symmetry vector. 
Aq, Bq and Ai, Bi(l — 1, 2, • ■ • , N) are used to denote the sites of carbon atoms. 
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Figure 2. An infinite atom chain and an atom Q out of the chain. Many atoms 
distribute regularly in a line form the infinite atom chain. The interval between 
neighbor atoms in the chain is T, and the site of atom Q relative to atom can 
be represented by numbers C\ and C2. 




Figure 3. The potentials inside the (20, 20) nanotube calculated from Eq.Q (squires) 
and the fit curve (solid curve). 
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Figure 4. The specific heats per atom c v of Ne inside (20, 20) nanotube at different 
temperatures without considering the 6, z effect. The triangles is the results of OS and 
the solid line is that of CSM. 
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Figure 5. The specific heats (triangles) per atom c v of Ne inside (20, 20) nanotubes 
at low temperatures. 
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Figure 6. A potential and its characteristic temperature. U m is the minimum of the 
potential and T c is the characteristic temperature. The parameters R\ and i?2 satisfy 
U(R 1 ) = U(R 2 ) = U m + T c . 
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Figure 7. The specific heats per atom c v of Ne inside (20, 20) nanotube at high 
temperatures. The triangles is the results of OS without considering the 8, z effect and 
the solid line is that of CSM with the 8, z effect. 



